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Abstract 



We review two novel techniques used to calculate tree-level scattering amplitudes effi- 
ciently: MHV diagrams, and on-shell recursion relations. For the MHV diagrams, we 
consider applications to tree-level amplitudes and focus in particular on the JV = 4 
supersymmetric formulation. We also briefly describe the derivation of loop ampli- 
tudes using MHV diagrams. For the recursion relations, after presenting their general 
proof, we discuss several applications to massless theories with and without super- 
symmetry, to theories with massive particles, and to graviton amplitudes in General 
Relativity. This article is an invited review for a special issue of Journal of Physics 
A devoted to Scattering Amplitudes in Gauge Theories. 
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1 Introduction 



There has been rapid progress since 2004 in the study of new structures in scatter- 
ing amphtudes in diverse quantum field theories, and their relationships with string 
theory and twistor theory. This has frequently come about via sudden and surprising 
insights into new properties, or new understandings of known properties, of tree am- 
plitudes - a subject that one might have been expected to have been fully understood. 

The initial work, following Witten's pioneering investigations into the twistor 
space localisation of amplitudes studied further the localisation of tree ampli- 
tudes. An early and vital insight was that one could glue together MHV tree ampli- 
tudes in order to obtain non-MHV amplitudes this expansion of tree amplitudes 
in terms of MHV diagrams provided a very natural way of understanding localisation 
on collections of lines in twistor space for more general tree amplitudes^ This "MHV 
diagram" approach was then applied successfully to one-loop amplitudes in A/^ = 4 [8], 
J\f = 2, Af = 1 [9|[T0] and even pure Yang-Mills theories [llj. Important steps towards 
the understanding of the MHV diagram approach arose from applying a particular 
non-local field redefinition to the lightcone Yang-Mills action [T2l[T3]. Extensions to 
non-MHV loop amplitudes have proved more difficult and this has impeded progress 
in this direction. More recent work has found very interesting and related new struc- 
tures in the integrands of loop amplitudes, by reformulating them in momentum 
twistor space or dual momentum space. Some of these topics are described elsewhere 
in this volume. 

A second, and particularly striking new discovery was that of recursion relations 
for gauge theory at tree level [TUITS]. These followed from a beautiful and elegant 
analysis of the complex analytic structure of tree amplitudes, together with the basic 
physical principle of unitarity of the S- matrix [16] . It was surprising that these recur- 
sion relations had not been found previously; they epitomise a concrete realisation of 
the old S-matrix programme, by providing a method to compute S-matrix elements 
directly from the analytic properties of the amplitudes, without ever referring to a 
Lagrangian, and requiring only on-shell quantities as input. These recursion relations 
have proved to be a very powerful and efficient tool as one may build more complex 
amplitudes from simpler ones. They have also been successfully extended to cover a 
wide class of field theories, including gravity [TTtfTS] and also loop amplitudes [19H2I], 
enabling new results to be derived and new structures to be found. Since they follow 
from basic properties of quantum field theory, interesting new applications continue 

^Note that in the context of twistor string theory the MHV diagrams are naturaUy related to 
the so-caUed disconnected prescription, while Witten's original approach [1 , further studied and 
extended in subsequent papers [3H5], was termed the connected prescription. Both prescriptions 
have been shown to be equivalent in [6l, where also intermediate prescriptions were identified, see 
also [7]. 
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to be found - recent examples include recursion relations for loop integrands in A/" = 4 
super Yang-Mills (SYM) [22||23] and recursion relations for partially off-shell quan- 
tities such as form factors [2l] - and it is to be expected that these techniques will 
hereon be an essential part of the armoury available for tackling problems in quantum 
field theories. 

The present chapter of this review is organised as follows: In the next section we 
will review the approach to tree-level amplitudes based on MHV diagrams, describing 
the method with some simple examples, as well as describing how it can be applied at 
one loop. In Section 3 we turn to recursion relations, with examples of how they can 
be used in gauge theories with massless and massive particles and in gravity. Finally 
we describe the manifestly supersymmetric formulation in A/" = 4 SYM using Nair's 
on-shell superspace [25j . 



2 The MHV diagram method 



The key quantity in the MHV diagram method is the MHV scattering amplitude of 
n gluons, given by [26ti29] 



/Imhv(1+, 2+, . . . , r , . . . , . . . , n+) = z(27r) V"'5^'^ ( ^^^0 ^'^^ 



(12)(23)---(nl) 

(2.1) 

A remarkable feature of this strikingly compact expression is that it is a holomorphic 
function of the spinor variables of the n scattered gluons. In QJ, Witten exploited 
this fact to perform a Fourier transform of this amplitude to Penrose's twistor space, 

Yld%e'^^^'^^^AMKW = Apt d'^x Hd^^^i^t + ^'"Ka) , (2-2) 
i=i ^ i=i 



where ^ 

■= (12)(23)...(nl) • ^^-^^ 

The right-hand side of (12. 2p vanishes unless each twistor coordinate of the n gluons 
Xa,i, fJ'a,i satisfies the relation 

fx^ + x^^Xa = 0, d = 1,2 . (2.4) 

This is the equation of a line in twistor space; in twistor theory, this is a famous 
relation - the incidence relation, which associates points in conformally compactified, 
complexified Minkowski space to complex lines, or CP^'sin twistor space [30] . Keeping 
this observation in mind, Cachazo, Svrcek and Witten (CSW) suggested in p] that 
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one may think of an MHV amplitude as an effective, local interaction in spacetimejj 
The idea is to introduce a particular off-shell continuation of an MHV amplitude, the 
MHV vertex. By gluing these MHV vertices in an appropriate way, to be discussed 
momentarily, we will construct amplitudes with an arbitrary number q of negative 
helicities - what we will call an N'^~^MHV amplitude. 

We can easily count the number of vertices v that are needed in order to calculate 
N''~^MHV amplitude, v MHV vertices provide us with 2v negative helicities, since 
each vertex has the MHV helicity configuration. At / loops, a generic diagram has 
V — 1 + I internal lines, each of which absorbs exactly one negative helicity. This is 
because a propagator connects an outgoing gluon of a certain helicity on one MHV 
vertex with an outgoing gluon of the opposite helicity on another MHV vertex. We 
are thus left with q = v + 1 — I negative helicities, and hence 

v = q-l + l_ (2.5) 

At tree level (/ = 0) we will then draw all possible MHV diagrams with q — 1 MHV 
vertices and "dress" them with the external gluons in such a way that each vertex 
has the MHV helicity configuration. 

Now two important issues have to be addressed in order to build up the diagrams. 
Firstly, vertices must be connected with an appropriate propagator. Secondly, (12.11) 
gives the expression for an on-shell amplitude, and, therefore, one has to provide an 
off-shell continuation of such an amplitude in order to use it as a vertex. In [2] , CSW 
provided us with answers to both questions, as we now describe. 



2.1 Propagators and off-shell continuation 

The answer to the first question just posed is extremely simple - the internal propaga- 
tors connecting MHV vertices are simply scalar propagators, i/iL"^ + is). This result 
can be justified from the Lagrangian formulation of MHV rules described in |12j . 
Next, we consider the issue of off-shell continuation. Let L be an internal (hence 
off-shell) momentum. The MHV amplitude depends on the holomorphic spinors as- 
sociated to the external, on-shell particles, and we will now provide a prescription 
that associates a spinor to the internal leg with (off-shell) momentum L. To this end, 
one introduces an arbitrary lightlike reference vector rjaa = VaVa- Then, there is a 
unique way to decompose L as [HET] 

L = I + ZT] , (2.6) 

^Mansfield showed in [12] that the MHV rules, to be explained below, can be derived at tree 
level by a canonical change of variables in the Yang-Mills path integral. The resulting action is local 
(only) in lightcone time. We briefly review this approach in Section 
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where / := Xl^l is a null momentum and 2; is a real parameter, which can easily be 
seen to be equal to ^ = L^/{2L-ri). The off-shell continuation proposed in [2] consists 
in using the spinor Al as the off-shell continuation for the internal leg of momentum 
L. It then follows that 

^L,a = (2.7) 

These equations coincide with the prescription originally proposed in [2] by CSW for 
determining the spinor variable associated with the off-shell (i.e. non-null) four- 
vector L defined in f l2.6p . The denominators arising from the right-hand side of fl2.7p 
will be irrelevant for our applications, since each MHV diagram is invariant under 
rescalings of the internal spinor variables; hence we will discard them and simply 
replace XL,a LadV°'- Each MHV diagram will therefore depend on the particular 
choice of fj; one will then have to prove that summing over all MHV diagrams, this 
spurious dependence cancels, showing that the MHV diagrams produce covariant 
answers. In [2] this was originally shown for diagrams involving one propagator 
i.e. NMHV gluon amplitudes. Later an elegant proof was found in [32] where it 
was shown that MHV diagrams are in direct correspondence with on-shell recursion 
relations (which by construction lead to covariant answers and will be described in 
later sections) where all negative helicity gluons are shifted. This was later extended 
to include all tree amplitudes in A/" = 4 S YM [33] . 



2.2 Examples of application 

The proposal of [2] that tree amplitudes could be built by gluing together MHV tree 
amplitudes as vertices, joining them with propagators, was very intuitive and gave an 
answer to the basic question of how tree amplitudes for n-particle scattering could be 
relatively simple, when they were obtained by summing a rapidly increasing number 
of Feynman diagrams. This method promised a new, efficient means for calculating 
amplitudes, seemingly sidestepping the complications of gauge fixing and dealing with 
unwieldy sets of diagrams. 

We will start with some simple applications at tree level and defer the more 
complicated example of all NMHV tree amplitudes for later where we will combine 
the CSW rules with Nair's on-shell superspace formalism. 

As a warm-up we will now rederive the well-known fact that the four-point tree 
amplitude (1"'"2~3~4~) vanishes. More generally, all trees with equal helicity gluons 
or only a single gluon of opposite helicity are zero. A very important exception to 
this are the three-point amplitudes (1+2^3^) and (1^2+3+) which vanish for real 
momenta in Minkowski signature but are non- vanishing for complex momenta (as we 
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will see later, this has profound consequences). Formally, our example is a so-called 
next-to-MHV (NMHV) amplitude, as it has three negative helicity gluons. 

There are two diagrams contributing to this amplitude, where two three-point 
MHV vertices are connected by a propagator, see Figure [U Using the fact that 




Figure 1: The two MHV diagrams contributing to (1+2 3 4 ) 

p2 = +P2y = (12) [21] and P'^ = {pi + Pif = (14) [41] we obtain for the sum of 
the two diagrams, 

.3 {2P)^ 1 (34)3 .3 (23)3 1 {P'A)^ 

' (PI) (12) (12) [21] (4 -P){-P3) (3 -P'){-P'2) (14) [41] (41) (IP') ' ^ ' 

Now we have to use the off-shell continuation introduced earlier namely Xp^a = 
PaaV'^ = —iX-p^a, and we have denoted incoming internal momenta with an ad- 
ditional minus sign since they correspond to an outgoing momentum —P. In our 
convention for on-shell momenta, Paa '■= AoAq, whenever we flip the orientation from 
outgoing to incoming we multiply A and A with i. Taking into account these extra 
factors of i and using the off-shell continuation we get, for the first diagram|§ 

(2|P|r/]3(34)3 _,(34) [lyf ^^.g) 



[7^|P|l)(12)2[21](4|P|r^][r^|P|3) [21] [27^] [37^] [4.^] 

where we used P = —pi—p2 = Ps+Pa and identities such as (2|P|?7] = {'2\—pi—p2\r]] 
— {21)[1t]] — (22) [27]] = (12)[1?7]. The evaluation of the second diagram gives 



.(14) [ly? 

'[23] [2r/][3r/][4r?]- ^ ' > 



^Note that, as we mentioned earlier, we are allowed to drop the normalisation constant in the 
definition of the off-shell continuation since ()2.9p is homogeneous in the internal ofF-shell spinors - 
this is a generic feature of MHV diagrams. 
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The sum of (12.91) and (I2.10p vanishes, since 



(34) (14) _ [23](34) + [21](14) _ [2|p3+Pi|4) _ [2^2+^414) _ 

[21] [23] [21] [23] [21] [23] [21] [23] ' ^ ' ' 

confirming that the (l+2~3~4~) amphtude is indeed zero. 

This procedure can be easily generahsed to calculate arbitrary n-point amplitudes 
with k negative helicity gluons. Such N'^'^MHV amplitudes require one to sum all 
possible MHV diagrams with k — 1 MHV vertices connected by scalar propagators. 
Applications are not restricted to amplitudes with gluons only but can also involve 
fermions or scalars [341135] . 



An obvious next step was to see if this method could be used to calculate loop 
amplitudes and we would like to briefly describe how this is achieved p]. Much was 
already known about loop amplitudes, particularly after work using unitarity methods 
in the 1990's |36lflT] . Applications of modern unitarity techniques are discussed in 
great detail in the chapters of this review by Britto ^2] and Bern and Huang 



The first and simplest example of a loop amplitude is the one-loop MHV amplitude 
in the M = 4 theory, which is given as sum over "two-mass easy box functions" 
Explicitly, one has 



n [fl-l 



(fc L 2 J 1 



i=l ■r=l 



In order to define the two-mass easy box functions i^^me^ ^^^^ introduce the scalar 
box integral I4, defined as 

h = -^(477)2 ^ j j—^^ p2(p _ K^Y^p j^^y^p + j^^y ■ (2-13) 

In the two-mass easy case, two opposite momenta in the box integral are massless 
{p and q in Figure |2] (b)), and the remaining two opposite momenta are massive 
(P = /cj H — ■ + ki+r-i and Q = /cj+r+i + ■ — h in the same figure). Denoting by 

/|.™f the two-mass easy box integrals, the corresponding F functions are defined by 



J-A- 



9 Tp2m e 

2me '^^n:r;i 



Jr+l],[r+l] _M [n-r-2\ 
''j ''i ''i+r+l 

Im 



2F" 

/I™ = IE— (2 14) 

where tf^ = (fcj + ■ ■ ■ + fej+^-i)^ are the kinematic invariants constructed from sums 
of cyclically adjacent external momenta. Note that = tf^ and = 
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(a) 



P 



(b) 



Figure 2: On the left we represent a two-mass easy box function F^^.^. n is the total 
number of external legs, the labels i and r are defined in the figure. On the right we 
depict the same box function in a slightly simplified notation, used later in (12.161) . 
The massless legs are called pi-i := p and pi+r '■= 



In the second equation in ( 12.141) we have considered separately the one-mass case, 
i.e. the special case where one of the two massive corners P oi Q becomes mass- 
less. The corresponding function is obtained as a smooth limit of the two-mass easy 
configuration [H]. Introducing the variables 



[P + pf = t. 



[r+l] 
i-1 ' 



t := (P + qY = t 



[r+l] 



(2.15) 



the explicit expression of the box function is 



F{s,t,P',Q' 



+ L^l--)+Ll.(l-- 



Llo 1 



{{-Qf) 



Ll2(l 



- L12 1 



St 



t 

(2.16) 



where P + p + Q + q = ^. The relation to the functions F^^.f is obtained by setting 
p = pi_i, q = pi+r, and P = Pi^ h Pi+r-i- 

Following the CSW proposal discussed above, one would immediately expect that 
the one-loop MHV amplitude should arise from the one-loop MHV diagrams depicted 
in Figure [3] below, i.e. from gluing two MHV vertices together with two propagators in 
order to form a one-loop diagram, and summing all possible diagrams while preserving 
the cyclic ordering of the external legs. This calculation was performed in [8|. One 
important point to discuss is the loop Integral measure. 



dM 



d^Lo 



L\ + is Ll + is 



SW^L2-L, + Pl 



(2.17) 
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Figure 3: One-loop MHV Feynman diagram, using MHV amplitudes as interaction 
vertices, with the CSW off-shell prescription. 



where Li and L2 are loop momenta, and Pl is the external momentum flowing outside 
the loopj so that L2 — Li + Pl = 0. Note that we also include the two propagators 
in the definition of the integration measure for convenience. 

In order to calculate the one- loop MHV diagram in Figure [31 we need to re-express 
a standard loop momentum in terms of the new variables / and z introduced in fl2.6p . 
One finds that [8] 

^ = d^l) ^ , (2.18) 

with the Nair measure [251 



dATil) := (/ dl)dH - [ldl]dH . (2.19) 

Note that the product of the measure factor with a scalar propagator d'^L/L'^ of 
(I2.18P is independent of the reference vector rj. Also, the Lorentz invariant phase 
space measure for a massless particle can be expressed precisely in terms of the Nair 
measure: 

dn5^^\f) = (2.20) 

where, as before, we write the null vector / as laa = lala, and in Minkowski space we 
identify / = ±/* depending on the sign of the energy. 

We are now in position to write the one-loop integration measure in terms of 
variables that are appropriate for the MHV diagram method. Expressing Li and L2 
as in (1^ . 

Li-^a,a ha^ice "I" ^iVaVa y ^ I5 2 , (2.21) 

we rewrite the argument of the delta function in fl2.17p as 

L2-Li+PL = l2-h + Pl;z , (2.22) 



^In our conventions, all external momenta are outgoing. 
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where Pl-^z '■= Pl — zr], and z := zi — Z2- Notice that we use the same rj for both 
the momenta Li and L2. Using fl2.2ip . we can re-cast fl2.17p as |8] 



dzi dz2 

dM 



Zi + iSi Z2 + 162 



d% d% 
2/10 2Z20 



5^'\l2 - h + Pl; 



(2.23) 



where := sgn(r7o/io)^ = sgn(/jo)£^, i = 1,2 (the last equahty holds since we are 
assuming r^o > 0). 

The next step involves integrating out the variable z' := zi + 2^2, which only 
appears through the integration measure 

dzi dz2 dzdz' 

(2.24) 



[zi + isi) {z2 + 162) ' [{z' + z + iei){z' - z + iea)] 

This can be done using the residue theorem, with the result [811^ that one can simply 
replace 

dzi dz2 ■ dz 
— — 2m -. (2.25) 

[Zi + lEij [Z2 + 162) Z + ie 

The last crucial step consists in performing the integration over z. Perfectly in tune 
with the S-matrix programme, this was converted in [8J into a dispersion integral. 
The final result for the integration is then [51 

dM = 2711 e{PlJ -2 gl dUPS{q,-l^;PL;z) , (2.26) 

where 

tiLIPS(/2",-/+;Pi;,) := d%6^+\ll) d%6'^-\ll) 6^^\12-Ii + Pl,z) (2.27) 
is the two-particle Lorentz invariant phase space (LIPS) measure, and we recall that 

5^1^) := e{±iMf). 

It is important to observe that the integration is performed for P|.^ > 0. By 
setting all the various external kinematical invariants P| to negative values, no poles 
are encountered along the integration contour and the is prescription can be dropped. 
However, f l2.26p provides us with the correct analytic continuation to the physical 
region, which is obtained by simply performing the substitution P| — > P| + ie. 

Having fully re-expressed the integration in terms of MHV variables, one can 
now perform all the remaining integrations and sum over all possible MHV diagrams 
corresponding to the different choices of particles running in the internal legs on each 
of the two MHV vertices. The interested reader can consult |8] for the remaining part 
of the calculation; in the following we will outline its main steps. 
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Firstly, one finds that each MHV diagram leads to a sum of four terms which are 
written as dispersion integrals in the channel fixed by the diagram. These are related 
to unitarity cuts of four different box functions. Crucially, in order to reconstruct 
complete box functions one has to take an appropriate combination of four dispersion 
integrals from four different MHV diagrams (channels). Surprisingly, it was found 
in [H] that dependence on the reference spinor 77 cancels out between the four dis- 
persion integrals contributing to a single box function and not only in the sum over 
all diagrams. The covariance at one-loop for arbitrary amplitudes was later shown 
in [l5] using the Feynman tree theorem. Collecting all the box functions, one finds 
precisely the expected result fl2.12p . 

We would like to make a few comments before closing this subsection. 

Firstly, the CSW approach, and hence the general MHV diagram proposal for 
calculating scattering amplitudes, was originally motivated by the twistor space lo- 
calisation of amplitudes and Witten's formulation of A/" = 4 super Yang- Mills theory 
as a theory in twistor space. At first sight, then, one might not expect these methods 
to work for theories with less supersymmetry (at tree level of course the gluon scatter- 
ing amplitudes are the same for the supersymmetric and non-supersymmetric cases). 
Nevertheless, one could just put aside considerations of the twistor picture, and sim- 
ply see if the MHV diagrams do yield the correct amplitudes in Yang-Mills theories 
with less than maximal supersymmetry. This was done in [9l[T0], where results from 
the MHV diagram calculation in A/" = 2 and M = 1 super Yang-Mills were found 
in agreement with those derived in [37] using unitarity. This success was somewhat 
unexpected, and led to a new, simpler understanding of the twistor space localisa- 
tion of loop amplitudes in these theories [IHlllZl- We also point out that this success 
was not accidental. Indeed, in |3S] it was proved that any non-MHV amplitude in 
supersymmetric gauge theories is reproduced correctly by MHV diagrams. 

Finally, the non-supersymmetric Yang-Mills theory was studied at one loop using 
MHV methods [TT]. Again, this correctly yielded the (cut-constructible part of the) 
amplitude, and this work was notable in that it included the first calculation of a 
previously unknown amplitude using these methods. It was at this point that the 
issue of so-called rational terms arose - these are terms in amplitudes that are not 
generated by unitarity cuts performed in four dimensions. This problem is not present 
in supersymmetric Yang-Mills theories, but does appear in the non-supersymmetric 
theory. The MHV diagram approach, as used in [TT], reproduced correctly all the 
non-rational terms, but not the rational ones. Note however that the MHV vertices 
that one inserts into the diagrams are essentially four- dimensional objects, although 
one does use dimensional regularisation in the loop integral; it is expected that it 
is this fact that underlies the missing rational pieces, although this remains to be 
investigated fully. Calculation of the rational terms has been achieved by various 
methods that are detailed e.g. in the chapter of this review from Britto [42] . 
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2.3 Super symmetric MHV diagrams 



It is generally the case that formalisms where symmetries are explicit are more pow- 
erful both conceptually and for calculations. It is easy to generalise MHV diagrams in 
A/" = 4 SYM to a manifestly supersymmetric version. This is particularly convenient 
when one is interested in amplitudes with external fermions or scalars. The on-shell 
superamplitude formalism introduced in the following is also reviewed in the chapter 
by Elvang, Freedman and Kiermaier [38]. 

The J\f = A theory contains two gluons G'^lp) with helicities 1,-1, four fermions 
ipA of helicity +1/2, transforming in the fundamental of the R-symmetry group 
SU{4:)r of the J\f = 4 theory, four Weyl fermions with helicity —1/2 in the 
anti-fundamental representation, and six real scalar fields (corresponding to particles 
of zero helicity) 0[ab] in the antisymmetric tensor representation of the R-symmetry 
group. Here A, B = 1, ... ,4 are fundamental SU (4)r indices. 

Following Nair, one introduces four auxiliary Grassmann variable^ r^^, and com- 
bines all the physical states of the theory into a super-wavefunction, 

A B ABC 

$(r^,p) := G^{p)+vUA{p) + ^<l>[AB]{p) + eABCD'^^-y^ 

(2.28) 

All amplitudes with a fixed total helicity are then collected into a single object called 
a superamplitude. The supersymmetry charges q^, qa satisfy the algebra {q^, qBa} = 
Paa^B where, on shell, Paa '■= ^a^a- In terms of the r^'s these are realised as 

n ^ d 

(la = ^-^^"^ ' (lA- = ~^'^s-A ■ (2-29) 

1=1 i=l 

One can realise the g-supersymmetry manifestly on the amplitudes by pulling out a S- 
function of supermomentum conservation, 5^^\Yl'i=i ^iVi)- The MHV superamplitude 
is then expressed as |25] 

n 71 n _ 

A,MHv(l,...,^) := ^5^'nE^^^*)^^'HE^^^^)n7^^^ ' (2.30) 

with A„_|_i = Ai. As one can see from (12.281) . in order to scatter a particular external 
state with helicity /ij, we need to expand (I2.30p and pick the term containing pi = 
2 — 2hi powers of rji. For example, the gluon MHV amplitude (12. ip with negative 
helicity gluons i~ and j~ arises as the coefficient of rifrij in the expansion of (I2.30p . 

As an example, we will now apply this formalism to derive a compact formula for 
all NMHV tree-level amplitudes in A/" = 4 super Yang-Mills. In this case we have 



'Not to be confused with the reference spmor of the same name introduced in Section [53] 
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to connect two MHV superamplitudes with an appropriate off-shell continuation of 
the internal leg P. The off-shell continuation for the spinor variables is as usual 
Xp-a = Paa^"' = —iX-p-a whcrc ^ is the reference spinor. Both MHV supervertices 
have degree 8 in 77 while the NMHV superamplitude has degree 12. This suggests 
that we have to augment the super-MHV rules with an jd'^rjp integration for every 
internal leg P where rjp = —irj^p are the fermionic variables for the internal legs. It 
is easy to see that in terms of component amplitudes this guarantees that the correct 
helicity state propagates between the two MHV vertices while it reduces the total 
Grassmann degree to 16 — 4 = 12 as expected. 

We now move on to the calculation. There is only one type of diagram, depicted 
in Figure m below. Note that we do not assign specific helicities to external legs as we 
are considering superamplitudes and that for a superamplitude only the total helicity 
Yli=i fixed. Evaluating this, we are directly led to 




Figure 4: The super-MHV diagrams contributing to the NMHV superamplitude. 

A,NMHV = j ^^'7p-4.MHv(^ + 1, • • • , j)-^v4.MHv(j + 1, . . . , i) 

(zz + l)(jj + l) 



„-3 



1 



1 



nLi(^^ + 1) {3P){p^ + -p){-P3 + 1) p' 



X J^,(2.31) 



where the sum is taken over all inequivalent MHV superdiagrams, and 



left J 



(2.32) 



A=l k=i+l 
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with gieft = Yl=i+1 ^kVk, bright = Z)Li+l ^kVk, Qtot = gieft + bright, and P = fcj+iH h 

ki. Note that in the last hne we factored out an overall supermomentum conservation 
delta function and performed the t]p integrations. Hence, we find 

where we have factored out the n-point MHV super amplitude (I2.30p . 



2.4 Lagrangian derivation of MHV diagrams 



We have discussed above how MHV diagrams correctly yield amplitudes in Yang- 
Mills theories. If one prefers the view of the "constructive S-matrix" approach, which 
uses primarily analytic properties and on-shell quantities as input to calculate the 
scattering matrix, then one might conclude that the methods discussed above, and 
in the following section, are perfectly adequate. In fact it turns out that they are 
superior in many ways to traditional approaches. Yet it is natural to ask if one 
can derive the MHV rules directly from a Lagrangian - essentially as some sort of 
Feynman diagrams. This might also have the advantage of illuminating how the 
quantum theory of MHV diagrams is best formulated in general, and also assist with 
exploring possibly interesting further structures and different kinds of perturbative 
expansions of the S-matrix. 

Such a Lagrangian description has indeed been formulated, as we will now de- 
scribe. The key steps were taken in [l2l|T3], and subsequent work can be found 
in |19]453] . A more recent article |5l] may also be consulted for further references; in 
the following summary we will use the notation of this paper. 

The starting point is to quantise the Yang- Mills theory in the lightcone gauge, and 
investigate whether a suitable canonical transformation can take one from the light- 
cone Yang-Mills Lagrangian to one whose terms precisely lead to the MHV vertices. 
We observe that such a MHV Lagrangian must therefore have an infinite number of 
terms. If one chooses the Minkowski space lightcone coordinates 

x = ^(t-x^), x = ^(t + x^), z = ^(x^+ix^), z = ^(x^ -ix^), (2M) 
V2 y/2 V2 V2 

with gauge-fixing condition A = 0, then, after elimination of unphysical degrees of 
freedom, the Yang-Mills action can be written in terms of positive- and negative- 
helicity fields A = and A = Az as 

S=^jdx y c/^'x (£-+ + £-++ + £—+ + £'—++), (2.35) 



13 



with 



£-+ = TiA dd-dd]A 



= -Tr (dd-^A) [A, OA] , 
+ = -Tr [A, OA] (dd-^A) , 
£'--++ = -Ti[A,dA]d-^[A,dA], (2.36) 

with S is a constant-x quantisation surface and c/^x = dx dz dz. 

The first two terms above, C^^ + £^++, describe self-dual Yang-Mills theory [55] . 
This is a free theory at tree level, and the only non-vanishing amplitude is at one 
loop, when all particles have positive helicity. Classically, therefore, one would expect 
there to be a canonical change of variables which transforms these two terms into a 
free Lagrangian. Formally, we seek a new field a (non-local) functional of A on 
the surface of constant x, such that -|- C can be written as a free theory, 

C-+[A,A]+ C-+^[A,A] = C-+[B,B], (2.37) 

with the condition on B that the transformation be canonical: 

**ae°(i.f)= / d'v i'^}!':^}M ''(i,y)- (2.38) 
After Fourier transforming, one can show that i3 is a power series in A, 



B{x,p) = Aix,p) + J2 / 

n=2 



(27r)3 ■ ■ ■ (27r)3 



n 

X 



-A[x, ki) . . . A[x, kn) 



(p, ki) (p, /ci /C2) • • • (p, ^1 H h A;„_i) 



(2.39) 



where we have introduced (i, j) := ij — ji = p\pi — P+pI- Furthermore, we note that 
^ is a power series in B, which is also linear in B, 

A{x,p} = B{x,p} + 

00 -m . ,3, ,3, I 



m=3 s=2 



(27r)3 ■ ■ ■ (27r)3 p 

(27r)3 6%p- J2 kr)B{x, h)... B{x, K) ■ ■ ■ B{x, k^) , (2.40) 

for certain known coefficients H'*^^(p, — fci, . . . , —km)- One can then see that when 
the lightcone action is written in terms of the new variables {B,B), there is the 
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free, quadratic term plus an infinite sequence of interaction terms, each of which has 
precisely two fields B, together with increasing numbers of fields B, 

C[A,A] = C-+[B,B] + C—+[B,B]+C—++[B,B] + C—+^+[B,B] + -- - . (2.41) 

Explicit calculation [19] shows that these interaction terms precisely give the Parke- 
Taylor amplitudes, continued off-shell in exactly the same way as proposed by [2]. 
Note that the reference null-momentum rj"'"' = i]"-!]"- thus comes to be identified with 
the gauge fixing vector in the lightcone action. 

We thus seem to have a candidate Lagrangian for generating the MHV rules. 
However, some care is needed here - to start with, the canonical transformation 
that we have used to change variables is non-local. This may affect the S-matrix 
equivalence theorem, by which we would like to conclude that the {B, B) variables 
may equivalently be used to describe the quantum theory. This theorem may in 
principle be violated by 1/p^ terms arising from the non-locality - in short, these 
could cancel LSZ reduction factors of and hence survive in the on-shell limit. 
In addition, and partly related to this, the Lagrangian above clearly cannot directly 
generate all amplitudes - e.g. the non-zero one-loop all-plus helicity amplitudes clearly 
cannot be obtained from the MHV rules. 

A more careful analysis reveals the following. Firstly, it has been found that the 
only diagrams that violate the equivalence theorem are certain dressed tadpole di- 
agrams [51] that explain e.g. the presence of the rational all-plus amplitudes, and 
apart from these one may freely use the {B, B) variables to calculate amplitudes. As 
far as the rational parts of other amplitudes are concerned, note that the transfor- 
mation of variables above can be defined in Z) = 4 — 2e dimensions and the resulting 
MHV vertices are also D dimensional - this differs from the MHV diagram method 
described earlier where four-dimensional vertices are used. The regularisation em- 
ployed may be the four- dimensional helicity scheme [S5] , where the internal momenta 
are in D dimensions whilst the external particles carry four- dimensional helicities, as 
in [19], or standard dimensional regularisation as in |50]. In either case, the additional 
e-dependent terms arising from these vertices are expected to combine with divergent 
terms from loop integrations to produce the missing rational terms of the amplitudes, 
although this remains to be investigated fully. 

Another way to produce the rational terms has been described in [53]. This 
is based on the lightcone gauge formulation of Yang-Mills theory given in [57H59] . 
which uses a regularisation in four dimensions most naturally formulated in terms of 
region (or T-dual) momenta. In this approach, a very simple two-point counterterm 
reproduces all the n-point all-plus amplitudes. It appears plausible that all other 
rational parts of amplitudes will arise in a similar way, although this has not been 
shown explicitly. 
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In conclusion, a Lagrangian containing interaction terms that reproduce the MHV 
vertices has been shown to be derived from ordinary Yang-Mills theory. However, in 
general there is a rather subtle interplay of (violations of) the S-matrix equivalence 
theorem and the regularisation scheme used, which is responsible for generating the 
rational parts of amplitudes. Whilst the schemes mentioned above do correctly gen- 
erate these terms in the cases studied, the calculations involved appear to become 
difficult in full generality. 

We also mention that an alternative derivation of the MHV diagram method di- 
rectly in twistor space was found in where, building upon [511 E2], the MHV 
diagram expansion was derived from a particular axial gauge fixing of a supersym- 
metric twistor action. See also [63j for very recent related work. 



2.5 Other directions 

To conclude this section, we would like to mention some other important applications 
of, or contributions to MHV diagrams. 

1. Recent work [M] (see also j65|l66] ) has described a momentum twistor space [67] 
formulation of the MHV diagram approach which is, roughly speaking, dual to that 
described above, and in which dual superconformal quantities play a key role. In this 
approach, the MHV vertices are replaced by unity, whilst the internal propagators 
are represented by superconformal i?-functions well known from other work. In this 
approach, features such as symmetries are readily apparent, and one may relatively 
easily write down integrands for quite general amplitudes and study their properties. 
The relationship with Wilson loops in twistor space can also be explored. A discussion 
of the space-time analogues of this work has been given recently in |66j . 




2. In an interesting paper [68], Dixon, Glover and Khoze managed to apply the 

MHV diagram method to the effective Lagrangian describing the one-loop top quark 
contribution to Higgs plus multigluon scattering processes, in the heavy top limit. In 
this limit, such processes are conveniently described by a dimension five operator of 
the form H Tr {F^^F'^"), where H is the Higgs field. Instead of discussing directly the 
amplitudes arising from that operator, the authors of [HSI considered the interaction 
0TrFgj-) + 0TrF^gj-,, where is a complex field whose real part is equal to H, 
H := (j) + (j), and -F(a)sd stands for (anti) self-dual part of the field strength. Because 
of the last equation, the Higgs amplitudes can be recovered by summing those with 
one (j) and those with one (p. The key point is that the (p and amplitudes are 
amenable to a derivation through MHV rules, which makes their calculation much 
more efficient. This approach was extended to one loop for the same amplitudes 
in [69]. Further extensions of the MHV diagrammatic method include applications to 
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QED [To] and Yang-Mills with coloured massive scalars [7T|[72]. 



3 BCFW recursion relation 



As we have alluded to above, one of the main ideas in the analytic S-matrix pro- 
gramme is that of calculating scattering amplitudes from the knowledge of their ana- 
lyticity properties, possibly without even knowing the Lagrangian of the theory |16j . 
A second key idea is to study amplitudes as functions of complex variables. We note 
here in passing that also for twistor theory the most natural arena is complexified 
Minkowski space. The BCFW recursion relation incorporates neatly these two pre- 
cepts. It provides us with an algorithm to calculate efficiently, and in a recursive way, 
all tree-level scattering amplitudes for various theories satisfying certain prerequisites 
(to be reviewed shortly), based on the knowledge of the singularities of amplitudes. 
At tree level, these are simple poles in the two-particle and multi-particle kinematic 
invariants, associated to collinear and multiparticle singularities, respectively. Start- 
ing from the smallest building blocks, namely three-point amplitudes, all amplitudes 
can then be constructed recursively. This is a very powerful statement - for example 
this implies that one can reconstruct the S-matrix of General Relativity simply from 
the knowledge of the three-point graviton amplitudes. We remark that these are just 
the squares of the three-point Yang-Mills amplitudes, which is the basis for another 
very interesting line of enquiry, relating gravity amplitudes to Yang-Mills ones. This 
is reported elsewhere in this review. In order to discuss the recursion relations, we 
consider first the case of colour-ordered amplitudes in massless Yang-Mills theory, 
and later generalise to other situations. 

A key feature of the BCFW recursion relation is that of mapping a subset of simple 
poles of the tree amplitude into simple poles in a single auxiliary complex variable 
z, whose residues are then calculated from well-known factorisation properties of 
amplitudes. In its simplest incarnation, one proceeds as follows. 



3.1 Derivation of the recursion 

Consider a generic amplitude A{pi, . . . ,Pn), and select two legs for special treatment; 
without loss of generality we can choose these to be 1 and 2 (note that we can also 
shift non- adjacent particles but this would lead to recursion relations involving more 
terms). One then shifts the two momenta as 





(3.1) 
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The shifts performed in (13.11) are chosen in a particular form in order not to alter 
the momentum conservation condition. Furthermore, we would like to preserve the 
on-shell condition for particles 1 and 2, which is possible if pii] = p2T] = 0. In real 
Minkowski space there are no solutions to these constraints but in complex Minkowski 
space there are two solutions, r] = A1A2 and r] = A2A1, where Pi = AjAj, i = 1,2, as 
usual. Picking for example the first solution, we rewrite (13. ip as 

Pi{z) = Xi{Xi + ZX2) ■■= XiXi , p2{z) = {X2 - zXi)X2 := X2X2 . (3.2) 

We can define the complex function A{z) := A{pi,p2,P3, . . . ,Pn)- Since pl{z) = 
p\{z) = for all values of z, A{z) defines a one-parameter family of scattering ampli- 
tudes. Importantly, the shifts introduced earlier make the kinematic variables com- 
plex, hence we should regard A{z) as describing a scattering process in complexified 
Minkowski space. 

It is very easy to prove that A{z) is a rational function of z, which has only simple 
poles in this variable. Indeed, singularities in tree amplitude arise from a Feynman 
propagator which is going on shell. Let Pij = Pi + ■ ■ - pj be the momentum of a 
certain propagator. There are three possibilities: either leg one or two belong to Pij; 
or, either both legs, or none, belong to Pij. It is only in the first case that Pij depends 
on z, since in the other two remaining cases, this dependence is either not present 
or it cancels since pi + P2 = Pi + P2- Focusing on the first case, and assuming for 
definiteness that particle 1 belongs to Pij, we can write Pij = Pij — ZX1X2, and 



Zij 



p2 P.2 z — Z-- 



(3.3) 



where Zij is the solution of P^j = 0, namely 



P2 

^ . = — . (3.4) 

' (1|P.,-|2] ^ ^ 

All the poles in z of A{z) are given by (13. 4p . Importantly, we also know the residues 
at these poles. Indeed, since A{z) is a physical amplitude, the residue will be given 
by factorisation on the corresponding pole. Specifically, as z ^ Zij, we have 

A{z) ^ 5^ A^(z = z,,)^Al\z = z,,) , (3.5) 



hence 



ResA(;^)|_^^, =-J2Al{z = z.,)t^A,'\z = z,,), (3.6) 



p." 

h '3 



where we have used (13. 3p . In (13.50 and (13. 6p we have included a sum over all possible 
internal helicities, as required from factorisation theorems. 
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One can therefore consider the following integral over a large circle at infinity 




Here we have defined A^o '■= A{z — )■ oo). Clearly it is at this point that the particu- 
larity of each theory enters the scene. We can have theories such that one can always 
find shifts such that the boundary term A^o vanishes, and theories where this is not 
possiblell We defer the discussion of this important point for a moment, and for the 
time being assume that such a boundary term is absent. In this case, (13. 7p expresses 
^4(0), i.e. the amplitude we wish to calculate, as a sum of residues. Using (13. 6p and 
the condition Aoo = we recast (13. 7p as 



This is the BCFW recursion relation. It expresses the (unshifted) amplitude as a 
sum of products of two amplitudes evaluated at shifted, complex kinematics. The 
fact that two amplitudes occur in the recursion is inherited from the factorisation of 
the amplitudes, which itself is a consequence of unitarity. The sum is over all possible 
arrangements of particles such that only one of the two shifted momenta belongs to 
either the left-hand or right-hand amplitude. With colour-ordered amplitudes, this 
corresponds to a simple sum, as illustrated in Figure [5l 

Before concluding this section, we observe that the BCFW recursion relations had 
been originally found by a different route in [H]. In that approach, one starts from 
a particularly simple and interesting combination of infrared consistency conditions 
found in [7B]. This infrared equation expresses a tree-level amplitude in A/" = 4 
super Yang-Mills as a linear combination of the coefficients of the one-mass and two- 
mass hard coefficients in the expansion of the corresponding one-loop amplitude. 
These coefficients can then be written using generalised unitarity [77], and further 
elaborated using the explicit expressions of the three-point amplitudes appearing in 
their expression ^14j . 



3.2 Examples in Yang-Mills theory 

Our first example of an application of the BCFW recursion relation will be the cal- 
culation of the colour-ordered gluon amplitude y4(l~2+3^4+) in Yang-Mills theory. 

^Of course, in the presence of a pole at z — )■ oo one can still write down recursion relations as 
long as one knows what the boundary term A^o is. See references |73fl75j for work on this important 
issue. 




(3.8) 
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Figure 5: One of the recursive diagrams contributing to the BCFW recursion relation 
for a colour- ordered amplitude . . . ,n). The particles with shifted momenta are 
adjacent - namely 1 and 2. Therefore, there is a single sum in the recursion relation, 
labeled by j . 

The first point to discuss is the choice of the shifts. We can understand this from 
the Feynman diagram point of view. When we perform shifts we will produce a 
certain sequence of, say m, propagators, between the shifted particles which depend 
on z. Their arguments are of the form 1/(P + zrj)'^, which for large z behaves as 
[{P + zr])'^]~"^ ~ z~"^. Among all the diagrams, those which are leading at large z 
are those where only three-point vertices appear, each of which providing a single 
positive power of momentum - and hence of z at large z. There will be exactly m + 1 
such vertices, giving a factor of 2;™+^ at large z. Together with the propagators, this 
diagram behaves as z at large z. In order to obtain a scattering amplitude we have 
to contract external lines with polarisation vectors. For gluons these are of the form 

+ _ KVa _ >^aVa . . 

'aa ■ ^^^^ , ^aa ■ ^3^^^ • l^-^J 

The only polarisation vectors which depend on z are those of the shifted particles. It 
is clear that if we shift the holomorphic spinors A of a positive-helicity gluon, and the 
anti-holomorphic spinor A of a negative-helicity gluon, these will provide a further 
suppression of z^'^ at large z. In conclusion, the diagram with the worst behaviour 
at large z is suppressed as z~^ at large z, on the condition that we choose shifts as 
specified above. 

Therefore, we will chose our shifts as 

\i:=~Xi + z\2, \2:=\2-z\i, (3.10) 
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with Ai and A2 unshifted. In the hterature such a shift is often denoted as a [12) shift, 
and following the discussion from above we know that if legs 1 and 2 have helicities 
— h the shifted amplitude A{z) vanishes for large z as required. Using MHV diagrams, 

it can be shown [15] that the helicity choices ++ and lead to the same benign 

large- 2; behaviour. 



1" 2 




Figure 6: The two diagrams contributing to the Yang-Mills amplitude 
A(l~, 2+, 3^, 4+). The diagram with the internal helicity assignment (H — ) vanishes 
because of the shifts (13.101) . 

Now we move on to the calculation. There are two recursive diagrams to evaluate, 
depicted in Figure [6], differing in the helicity assignment of the internal gluon. In both 
diagrams, z is evaluated at the position of the pole (23) [23] = 0, or 

^ = ^. (3.11) 
(13) ^ ' 

Consider first the diagram with internal helicities (H — ). In this case the amplitude 
on the right-hand side has the MHV helicity configuration, and it is easy to see that 
it vanishes. Indeed, at the pole we have (23) = (3-P) = (P2) = 0. We now calculate 
the diagram with the opposite helicity assignment, ( — h). It is equal to 

J^\^(z!tM\ (3 12) 

(P4)(41) J \[23][3-P]y 

In our conventions, in which all particles have outgoing momenta, the spinors in the 
right-hand amplitude associated to the internal particle with on-shell momentum —P 
are I — P] = i\P]. We also have 

(1P)[P2] = (1|P|2] = (lb2+P3-^AiA2|2] = (13)[32] . (3.13) 
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In principle in f l3.13p we would have to insert the value (13. lip of z at the pole, but 
in this particular calculation the shift drops out since [22] = 0. Further multiplying 
and dividing by a factor of (1P)[P2] we form in the denominator the combinations 



(1P)[3P] = -(lb2 + P3-;2AiA2|3] = -(12)[23] , 
(P4)[P2] = -(4^2 +P3 - ^AiA2|2] = -(43)[32] . 



(3.14) 



Collecting terms, (I3.12p becomes 



(13)^[32]4 1 



(13)^ 



(3.15) 



— I 



(41) [23] (12) [23] (43) [32] (23) [32] 



= I 



(12) (23) (34) (41) ' 



which is the expected result. We remark that at no point in this derivation we 
needed to plug in the explicit value of z in (13. lip . This is due to the simplicity of 
this calculation; in general derivations one does need to insert the explicit value of z 
at the poles for each diagram. 



3.3 Extensions to other theories 

The derivation presented above is very general, therefore one might expect to be able 
to apply on-shell recursion relations in different theories. Specifically, at three stages 
we have made reference to a specific theory: 

1. We have considered colour-ordered amplitudes in Yang- Mills theory. 

2. We have considered massless particles. 

3. We have assumed that v4(z) — )■ as z — )■ oo. 

Let us analyse these three points in turn. Firstly, colour-ordering is just a technical 
point. We could have just as well applied the recursion to complete amplitudes, and 
nothing would have changed except that the sum over all recursive diagrams would 
have become a double sum. We have also considered shifts appropriate for massless 
particles, but one could well extend these to massive particles as shown in [781179]. 
For example, if particle 1 is massless and particle 2 is massive, one can consider the 
shifts 




(3.16) 



7/ 



aa 



mm) 



(3.17) 
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which has the properties r] ■ t] = pi ■ r] = p2 ■ r] = as for the massless case, which 
are needed to preserve the mass-shell conditions of the unshifted momenta. There is 
another possibility namely we could have shifted by 

ry'- = {p2~Xir~Xi . (3.18) 

As in the massless case, we now consider the deformed tree- level amplitude A{z) and 
use the fact that this function can have only simple poles in z coming from internal 
propagators i/P{zY going on-shell. This can only occur if particles 1 and 2 are on 
opposite sides of the propagator so that P = pj + . . . + pn + Pi and P{z) = P + zr]. 
The propagator becomes 

(3.19) 



P{zy — m? — m? + 2zP ■ rj 

Assuming that A{z) for large z and using Cauchy's theorem leads in a similar fashion 
as in the massless case to the recursion relation 

n 

A = Y.M^j)p2^2M^j) , (3.20) 

where Pj = pj + . . . + Pn + Pi, and Zj = {w? — Pj)/(2Pj - rj). Of course, we could have 
picked any other pair of external particles and even selected two massless particles 
with the usual BCFW shifts, or two massive particles - although this last option is 
rather inconvenient for practical calculations. 

As a simple illustration, we now consider the case of four-point amplitudes with 
two massive scalars. In order to do this, we need the expressions of the three-point 
amplitudes with two massive scalars and one gluon as seed amplitudes. These were 
determined from Feynman rules in [78], 



^3(P^*^^.^) = %^, ^3(p?,*-.pJ) = -^^, (3.21) 

{r]ik) [7]2k\ 

with Pi 2 = ""^^ and k"^ = 0. Here, 771^2 are two reference spinors needed for the 
definition of the polarisation vectors, and it can easily be seen that the massive three- 
point amplitudes f l3.2ip are independent of the choice of r]i 2- 

Let us now focus on A{pi , 2"*",p2 ), the amplitude of two positive helicity glu- 
ons and two massive scalars. We choose to shift the two gluon momenta ki^2 with 
rj = A2A", hence A2 and Ai are unshifted. There is only one recursive diagram to be 
considered involving two three-point amplitudes (13.211) with a massive scalar propa- 
gating between them. This yields 

A3(p+,l+,-P-) p, ' ^ A3(P+,2+,p2') . (3.22) 
P^ — 
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Choosing reference vectors r]i = k2 and ri2 = ki for the two three-point amphtudes 



one arrives at 



'(12)2((pi + A:i)2-m2) 



(2|pi|l](i|p2|2] 



and using Fierz identities and on-shell conditions one can simphfy this further to find 



Similarly the amplitude with two massive scalars and two gluons of opposite helicity 
can be obtained from an almost identical recursion relation 



Other tree amplitudes involving massive scalars have been calculated explicitly in 
[78l[79]. Such amplitudes have particular relevance for the calculation of one-loop 
amplitudes in QCD using unitarity. Due to the presence of rational terms, already 
mentioned at the end of Section [2^ it is not sufficient to consider unitarity cuts in four 
dimensions - these have to be performed in Z) = 4 — 2e dimensions, if dimensional 
regularisation is used. If the particle running in the loop is a scalar, then going 
to D dimensions effectively turns a massless scalar into a massive four- dimensional 
scalar whose mass squared /i^ has to be integrated over (/x^ is identified with the 
— 2e dimensional part of the loop momentum). Applications of this approach to 
calculations of one-loop QCD amplitudes are described in more detail in the chapter 
of this review from Britto [42] . 



3.4 Recursion relation in Einstein gravity 

It is very interesting to consider on-shell recursion relations in theories other than 
Yang-Mills, for instance Einstein gravity. Following from the BCFW work, this was 
first done in [T71IT8] . A general formula for the MH V amplitudes for graviton scattering 
was first given in [80] . Choosing particles 1 and 2 to carry negative helicity, while the 
remainder carry positive helicity, the result of [80] is 




(3.24) 




(3.25) 



g [12][n-2n-l] 1 
(In-l) N{n) 



n— 3 71—1 



M(l,2,3,-- - ,n) = z(12) 



n n (^^■)^ + ^(2'---' 



n-2) 



i=l j=i+2 



(3.26) 



where 



F 




(3.27) 
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and V indicates permutations. The recursion relations for gravity on the other hand 
yield the following solution for the MHV amplitudes [T7] . 



(12)6[H„_2] . . {2\l^ + ... + ls-l\ls 



+ P(^l,...,^„_2), (3.28) 



where 



2 (2«i)(2Z2)(«l«2)(«2«3)(«l«3) 

For n = 5 the product term is dropped from (13.281) . 

This is with the choice of reference legs 1~ and 2~. Subsequent work has expanded 
on these results, finding new and interesting formulae for the MHV amplitudes and 
proving equivalences between the different results [8TI - I86] . There are intriguing links 
to Yang- Mills theory - see for example [HZ] and references therein - and these are the 
subject of on-going research (see [HnilHHllHH] and references therein). 

We would like to illustrate how recursion relations can be applied to gravity in a 
simple example. We will focus here on the four-point MHV amplitude of gravitons 
which reads 

M(l-,2-,3+,4+) = (3.30) 
where N{n) := Y[i<i<j<n{'^ j) ■ ^^^^ consider the shifts 

Ai := Ai + 2;A2 , A2 := A2 - -zAi , (3.31) 

with A2 and Ai unshifted. It is immediate to see that, under these shifts, one has 
M(l~, 2^, 3'*', 4+) ~ as z ^ 00, hence there is no boundary term in the recursion 
relation. There are two diagrams to consider, one of which is depicted in Figure [71 
The other is obtained by swapping the labels 4 with 3. The three-point amplitude 
on the left-hand side of Figure [7] has the MHV helicity configuration, whereas the 
amplitude on the right-hand side has the MHV configuration. There is in principle 
another diagram where the left-hand amplitude is MHV and the right-hand one is 
MHV, but this is easily seen to vanish, since the right-hand MHV factor would be 
proportional to 

l3P] = m>=Bm = 0. (3.32) 
We now evaluate explicitly the diagram in Figure [71 It is equal to 

m4^)(1,2,3,4) = M(3)(P+,4+,i-)-^M(3)(-p-,2-,3+) , (3.33) 
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Figure 7: One of the two diagrams contributing to the recursion relation for the MHV 
amplitude M(l", 2^, 3+, 4+). The other is obtained from this by swapping legs 3 and 
4. 



where 

M(3)(P+, 4+,i-) = , (3.34) 

V [41][1P] 

M(3)(p-,2-,3+) = f^^-^ 

' (23)(3 -P) 

and P2 = (p^ +p^)2_ Using {i P) = {t\P\l]/[P 1], we get 

^e'a.2.M) = .(j|||M_. (3.35) 

The full amplitude is M^^) (1", 2", 3+, 4+) = M^ei (1, 2, 3, 4) + Mr^ej (1, 2, 4, 3). We 
conclude that the result of the recursion relation is 

It is easy to check that fl3.36p agrees with fl3.3Qp . thus reproducing correctly the 
expected four-point MHV amplitude. More examples of applications of the recursion 
relation are given in [iTlllH]. A thorough analysis of the large-z behaviour of gravity 
amplitudes is performed in 
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3.5 Manifestly supersymmetric recursion relations in = 4 
SYM 



It is very natural to formulate manifestly supersymmetric on-shell recursion relations, 
and this is straightforward to write down [92||93] . The interest in such generalisations 
stems from the fact that having an algorithm that calculates amplitudes while re- 
specting all symmetries of the theory is bound to produce expressions that are simple 
and compact. In particular, it was shown in [92] using the supersymmetric recursion 
relations described below that the tree-level S-matrix of N' = 4 SYM is covariant un- 
der dual superconformal symmetry and that this covariance is manifest term-by-term 
in the recursion relation. 

The main idea consists in supersymmetrising the momentum shifts of the standard 
BCFW recursion by accompanying them with corresponding supermomentum shifts. 
Consider again the shifts in (13. 2p . There we shifted Ai and A2. The sum of the 
supermomenta of the shifted particles, qi + q2 := ?7iAi -|- ?72A2 should also be invariant 
under the shifts, hence the shift of A2, A2 = A2 — zXi, must be compensated by a 
corresponding shift of r/i, namely 171 (z) = r^i + zri2. Therefore, the set of (super)shifts 

\i{z) := \i + z\2 , \2iz) = \2 - z\i , r]i{z) = r]i + ZT]2 , (3.37) 

preserves momentum and supermomentum conservations, as well as the on-shell con- 
ditions. This produces a one-parameter family of superamplitudes, 

Aiz) := ^({Ai, Ai, r/i}; {Ai, A2, Vih • • •) , (3-38) 

where the dots denote the unshifted momenta and supermomenta for the remaining 
n — 2 particles. The supersymmetric recursion relation follows from arguments similar 
to those which led to ( 13. 8p . We have 

^ = 5Z [d'VpMzp)^^R{zp) , (3.39) 
p 

where r^p is the Grassmann variable associated to the internal, on-shell leg with 
momentum P. The sum is over all possible P such that precisely one of the shifted 
momenta, say pi, is contained in P. The two superamplitudes are then evaluated at 
the solution zp of the equation P'^{z) = 0, where P{z) = P + 2A1A2. 

Notice that superamplitudes are characterised by the number of external particles 
and their total helicity, which is the sum of the helicities of all external particles. 
Hence, in the recursion relation (I3.39P we have an important constraint on Al and 
Ar, namely the total helicity of Al plus the total helicity of Ar must equal the 
total helicity of the full amplitude A. This condition replaces the sum over internal 
helicities in the standard BCFW recursion relation (13. Sp . 
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In writing f l3.39p we have assumed that A{z) — )■ as 2 — )• oo, as was demonstrated 
in [921193] . Here we will present the approach of [93], showing that 

Am=4{z) ~ - , Aj^=s{z) ~ ^ , as z -)■ 00 . (3.40) 

z 2;^ 

The key point is that one can use maximal supersjTumetry to set to zero two of the t]^ 
variables, say rji and 772, in a given superamplitude ^(Ai, Ai, rji, A2, A2, rj^, - ■ ■ ; A„, A„, ?7„) 
in maximally supersymmetric A/" = 4 super Yang-Mills or A/^ = 8 supergravity. In- 
deed, consider the particular combination of q supersymmetries defined by := Cs^f ; 
where B = 1, . . . ,Af and determine the 2A/' parameters by requiring that 

or, equivalently, q(ri^2 — ~Vi2- easily finds that 

= 1^ {-K<^V2 + k<^vt) , (3.42) 
hence the action on a generic coordinate rji is 



e^c^, := r,[ = r/, - 771 + r/2 1^ , (3.43) 



and in particular e^'^rji = e'^'^ri2 = 0. 



Invariance under q supersymmetry of an arbitrary superamplitude A, which we 
write as 

^:=5(4)(p)5(2A^)(g)A(Ai,Ai,r/i;A2,A2,r/2;--- ;A„,A„,r/„) , (3.44) 

implies that 6^^\p)6^'^-'^\q) [q^A] = 0, hence (omitting for brevity the 6^'^\p)S^'^^\q) 
factors) e^'^A = A. Acting with the q operator explicitly, we therefore find that 

A{\i, Ai, 0; A2, A2, 0; ■ ■ ■ ; Aj, Aj, i]'-; ■ ■ ■ ; A„, A„, 1]'^) 
= A{Xi, Ai, T]i; A2, A2, V2]-- - ; Aj, Ai, A„, A„, r]n) , (3.45) 

where t]'^ is given in (13. 43 p . 



The application of the above to prove the large-z behaviour of superamplitudes 
works as follows [93]. Consider the family of amplitudes A{z) obtained by performing 
the supershifts in fl3.37p . and use f l3.45p to set ?7i(z) and 772 (-2) to zero. Importantly, 
the necessary supersymmetry parameter turns out to be z- independent: 

= 7^ {-kaV2 + haVl) = 7^ {-k»V2 + k^vt) , (3-46) 



28 



where we used f l3.37p . This imphes that 



A{z) = ^(Ai, Ai,0; A2, A2, 0; ■ ■ ■ ; Ai, Ai,?7-; ■ ■ ■ ; A„, A„,?7^) , (3.47) 

where none of the Grassmann variables 77^' in fl3.47p contain z, i.e. the only 2;-dependence 

occurs through Ai and A2. As a consequence, the behaviour of A{z) at large z will be 
the same of that of a gluon amplitude in Yang-Mills, or graviton amplitude had we 
started with J\f = 8 supergravity, with gluons/gravitons 1 and 2 both with positive 
helicity (since we have set to zero the corresponding t] variables for particles 1 and 2). 
In the Yang-Mills case, such amplitudes are known to fall off as I/2; at large z jl5] . 
whereas in gravity the fall off is l/z"^ [91] j§ and ( 13.40p follows. 



3.6 Example: MHV superamplitude 

The simplest application of the supersymmetric recursion relation of [92l|93] is to the 
derivation of the MHV superamplitude 

^MHv(l, • • • ,nj - i{2Tr) — — — . (3.48j 

(12) ■■■(nl) 

Performing the supershifts in fl3.37p . we find that there is a single diagram contributing 
to the supersymmetric recursion relation, represented in Figure Ei In this diagram, 
the right-hand amplitude is always a three-point MHV superamplitude, whereas the 
amplitude on the left-hand side is MHV. The expression of the former was found 
in [9211931 and is 



A nO',\ Y9 ^4e(4). ^ ^ , ^W(r7i[23] +V3m 

-4mhv(1'2,3) = i{2n) 6^ > [pi + P2 + Ps) ^^^J [23] [31] ' ^ ^ ^ 

In [92] it was also explicitly shown that this three point amplitude is invariant un- 
der supersymmetry, as well as covariant under the dual superconformal symmetry 
introduced in [95]. 



We describe in detail the derivation of the four-point superamplitude; the gen- 
eralisation to higher numbers of points is entirely straightforward. In this case the 
amplitude on the left-hand side is a three-point MHV superamplitude, 

^MHv(l,2,3) = t{27i) 6^ >{pi+p2+P3) (12) (23) (31) ' ^ ^ ^ 



^Some of this benign large- 2; behaviour, which is not apparent from a Feynman diagram analysis, 
was observed already in [171 [HI [90]. See also [94] for a more recent investigation of the large- 2: 
behaviour of scattering amplitudes in various theories. 
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Figure 8: The single recursive diagram contributing to the recursion relation for the 
n-point MHV superamplitude. The amplitude on the right is anti-MHV, whereas that 
on the left is an {n — 1) -point MHV superamplitude. 



Notice that f l3.49p and fl3.50p can both be determined from symmetry considerations 
alone (up to an overall normalisation). For instance, ( I3.50p is found by requiring that 
it depends only on the holomorphic spinors Ai, A2, A3 and satisfies the three equations 

k ^MHv(l, 2, 3) = ^MHv(l, 2,3), i = 1, 2, 3 , (3.51) 
where the helicity operator for particle z is [I] 



After this short detour, we now apply the supersymmetric recursion relation f l3.39p . 
with 

A, = 5^Hp.+h + P) (3.53) 

(1P)(P4)(41) 

^ ^ S^'^Kp2+P3 - P) S^'^Kvpm + mjS -P] + V3[-P2]) 

[-P2][23][3 -P] 

Using the identity 

S^'\mX, + r],X, + r^pXp) 5^^\r^p[2'i] + r^^[^ -P] + r73[-P2]) 
= <5(«) ( f,:\^ 5(^)(^p[23] + -P] + %[-P2]) , (3.54) 
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as well as 

^Vik = ^Vi^i, ^Pi = ^Pi^ (3-55) 

i i i i 

we arrive at 

^(1,2,3,4) = ^<5(^)(5^p,)5(«)(^r7,A,)A(l,2,3,4), (3.56) 

iGL,-R i€L,R 

where 

PI3 (41)[23](1P)(P4)[P2][3P] y ^"^^ ^ ' ^ ^ ' ^ ^ ' 

It is easy to see that (IP) (P4) [P2] [3P] = (12)(34)[23]2, therefore we conclude that 

/l(l,2,3,4) = , , , ^, , , , . (3.58) 
^ ' ' ' ^ (12) (23) (34) (41) ^ ^ 

We have thus reproduced the expected supersymmetric MHV superamphtude (13.481) 
in the four-point case. The generalisation of the calculation showed above to an n- 
point MHV superamphtude is very simple. The only difference is that the amplitude 
on the left-hand side of Figure [8] will be an {n — 2)-point MHV superamphtude. The 
algebra is identical to that of the four-point example discussed above and leads to 
the expected result (I3.48p . 

We conclude with a few references to relevant other work. In [HS] , an explicit so- 
lution to the TV = 4 supersymmetric recursion relation was found - this is described 
elsewhere in this review. In TV = 8 supergravity it is also possible to present an explicit 
solution to the recursion relations, see [97]. Applications of recursion relations to ra- 
tional terms in one-loop QCD have been considered in [T9H2T] . and are also discussed 
elsewhere in this volume. The relation of recursion relations and twistor space has 
been considered in [9811102] . Finally, recursion relations in three-dimensional super- 
symmetric Yang-Mills and Chern-Simons matter theories have recently been written 
down in [T03| fT0i]. 
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